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ABSTRACT 



This paper examines the crushing characteristics of a web 
girder subjected to a local in-plane crushing load. First, a 
simplified model of the post-buckling deformation of the web 
girder is developed. Then this model is used to formulate an 
upper-bound solution for the plastic behavior of the deforming 
web girder. Finally small-scale crushing tests of extruded 
aluminum double-hollow forms are conducted to verify the 
deformation mode and the upper-bound solution. The web girders 
of interest for this research are the longitudinal web girders 
in a unidirectionally stiffened double hull (USDH) structure. 
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1 . 



INTRODUCTION 



The ability of ships to withstand damage- -whether from 
collision, grounding, or explosion--is an active concern of 
naval architects and ship owners. In 1985, the Hitachi Zosen 
Corporation of Japan presented the concept of an advanced double 
hull design which utilizes longitudinal girders connecting inner 
and outer hulls to form a cellular structure that wraps around 
the girth of a hull. 111 Among other advantages, this 
unidirectionally stiffened double hull (USDH) structure, shown 
in Figure 1, offers the potential of increased protection 
against local damage resulting from accidental loading. This 
potential has contributed to recent interest in the USDH 
structure for both commercial and naval ship designs. 121 

Damage to ships involves a complex interaction of 
structural members. While several authors have developed models 
of ship collisions, research concerning ship groundings is in 
the initial stages. Most ship damage models have dealt with 
conventionally framed ships which utilize a gridwork system of 
large transverse web frames intersected by smaller longitudinal 
stiffeners. In general, for both collisions and groundings, 
these transverse web frames are subjected to large compressive 
in-plane loads. 

McDermott et al . 131 developed a mathematical model of the 
various stages involved in a minor collision of longitudinally 
framed oil tankers. This model was applied to both single hull 
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and double hull tankers. McDermott assumed that the resistance 



offered by a buckled web frame remained constant during 
deformation of the web frame in the transverse direction. 
Chang 141 in his discussion of McDermott's paper recommended that 
the post-buckling resistance of web frames be further 
investigated because many structural members offer decreased 
resistance after an ultimate strength has been reached. 

Wierzbicki et al. [51 investigated the extent of damages 
involved in high energy groundings of tankers for a study 
conducted for the Joint M. I . T . -Industry Program on Safe Tankers. 
For this study, a simplified model was developed to predict the 
energy dissipated by the crushing of transverse web frames. 

This model utilized the assumption that the deformation of 
the transverse web frame and the hull plating which acts as a 
flange could be separated into two mechanisms . It was assumed 
that the web detaches from the flange and rolls into two 
cylinders of equal radii, one on either side of the point of 
load application. It was assumed that the hull plating deforms 
into a trough shape. A sketch of this model is shown in Figure 
2. Geometric compatibility of the connection between the web 
and the attached plating, however, was not maintained. 

Instead of simplifying the post-failure resistance of 
deforming webs by modelling the resistance as a constant, the 
plastic behavior of the web may be described by a non-linear 
spring in a global model of a ship structure. 
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For example, the USDH structure can be considered to be a 
sandwich structure with a definable crushing strength of the 
core, a finite shear strength in the direction of the 
longitudinal web girders, and negligible shear strength in the 
transverse direction normal to the web girders . This sandwich 
structure can be modelled as rigid-plastic membranes 
(representing the inner and outer hull plating) connected by 
vertical (crushing strength) and diagonal (shear strength) non- 
linear springs . This model is shown in Figure 3. 

The objective of this study is to develop the load- 
deformation characteristics for longitudinal web girders 
subjected to local in-plane crushing loads. These 
characteristics can then be used to describe the vertical 
springs of the model described above. 

The crushing of longitudinal web girders presents a problem 
similar to that of crushing tubular members . Several studies 
have been conducted which have explored the large deflection 
plastic behavior of thin walled tubular members. Reference [7] 
contains a review of several of these studies . Much of the 
analyses involving crushing of tubular members is focused on 
crushing due to axial loads. This class of problems lends 
itself to simple modelling due to the axi-symmetrical cross- 
section of tubular members . 

The crushing of longitudinal girders in the USDH structure, 
however, is similar to crushing of tubes subjected to lateral 
loads rather than axial loads. De Oliveira et . al. 1 ® 1 present 
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an analysis of crushing of circular tubes due to lateral loads. 
This class of problems is more difficult to analyze due to a 
lack of symmetry and due to the fact that a tube's cross- 
sectional shape changes during deformation. If the crushing 
load is localized rather than distributed, the distortion of the 
cross-sectional shape occurs locally and propagates. This 
propagation of damage is difficult to predict with out 
experimental results . 

The problem under consideration for this study resembles 
crushing of a multiple box column-type structure subjected to a 
localized lateral load. Although de Oliveira's work focused on 
crushing of circular tubes, the same approach of using simple 
methods of structural plasticity applied to a simplified 
deformation model is used for this study. 

First, models for two modes of deformation of a web girder- 
-local denting and global bending- -are developed based on some 
photographic evidence of the deformation of fully clamped plates 
that were subjected to in -plane impact loading. Next, upper- 
bound solutions for the plastic behavior of deforming web 
girders are formulated for each mode. Finally, theoretical 
predictions based on the upper-bound solutions are compared to 
the results from a limited number of small-scale crushing 
experiments . 
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2 . 



PROBLEM DEFINITION 



2.1. UPPER-BOUND THEOREM: An upper-bound approach will be 
used to find the crushing characteristics of a web girder 
subjected to a local in-plane crushing load. Calladine [91 
provides a concise statement of the Upper-bound Theorem: 

"If an estimate of the plastic collapse load of a body 
is made by equating internal rate of dissipation of 
energy to the rate at which external forces do work in 
any postulated mechanism of deformation of the body, 
the estimate will be either high, or correct." 

The postulated mechanisms of deformation to which Calladine 
refers to must be kinematically admissible. 

Two modes of deformation are postulated for this study. 
The first mode involves local indentation of the web girder. 
The second mode involves the development of a three hinge 
mechanism to induce global bending of the web girder between 
transverse bulkheads . These postulated deformation mechanisms 
are described in sections 3 and 4 . Assumptions concerning the 
geometry and boundary conditions of the deforming body, the 
material properties of the body and the applied external loads 
for this solution are described in the following sections. 
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2.2. STRUCTURAL UNIT : Figure 4 shows a section of the USDH 
structure consisting of rectangular cells . This cellular 
structure runs between transverse bulkheads and wraps around the 
girth of a hull. Due to symmetry it is sufficient to define a 
structural unit which consists of one longitudinal web girder 
with the attached inner and outer hull plating. This structural 
unit is shown in Figure 5. 

The depth of the web girder is "H" and the thickness of the 
web girder is "t H " . Because the inner and outer hull plating 
act as fully effective flanges, the breadth of the attached 
plating is equal to the girder spacing, "B". The thickness of 
the hull plating is "t p " . The length of the web girder is 
determined by the distance between transverse bulkheads . This 
length is much greater than the girder depth and spacing. The 
structural unit, therefore, can be assumed to be effectively 
infinite in the longitudinal direction. 

Due to symmetry, the attached hull plating flanges are 
constrained from in-plane movement in the X-direction along the 
longitudinal edges (X = ± B/2) . Both ends of the structural 
unit are considered to be fully clamped at the transverse 
bulkheads (Y = ± L/2) . 

2.3. LOADING CONDITION: The web girder is considered to be 
subjected to a local crushing load such as that delivered by a 
knife-edge indenter or a rigid punch normal to the longitudinal 
axis of symmetry. This loading is quasi-static, therefore 
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dynamic effects are ignored. The applied load has a magnitude 
of P B . This load causes the flange to deflect a distance A 
relative to the un-deformed region of the flange. The rate of 
deflection is denoted by A. 



2.4. MATERIAL PROPERTIES AND ASSUMPTIONS: The material is 

assumed to be rigid/plastic isotropic time-independent material 
which can be described by the flow stress <7 0 . Strain hardening 
effects may be included by adjusting the value of the flow 
stress based on the average plastic strains of the deforming 
regions during the deformation process . The average plastic 
strain may be obtained by integrating the strain over the length 
of the deforming region: 

L 

= -7 fe(z) dz . ( 1 ) 

0 

Bending stresses and membrane stresses are modelled as 
uncoupled. Deforming regions which experience high bending 
stresses are considered to have negligible membrane stresses. 
Likewise, regions which develop significant membrane stresses 
are assumed to have negligible bending stresses . 
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2.5. LIST OF NOMENCLATURE : This section provides a list of 

the symbols and parameters used in this paper. 

HULL PARAMETERS 

H: Depth of web; 

B: Longitudinal girder spacing; 

t„: Thickness of web; 

t p : Thickness of inner and outer hull plating 

L: Distance between transverse bulkheads; 

H*: H / t„; 

B* : B / t p ; 

t*: t p / t„; 

L*: L / H; 

MATERIAL PARAMETERS 

C7 0 : Flow stress; 

G ult : Ultimate tensile yield stress; 

e ult : Ultimate tensile strain; 



M c : 


Fully 


plastic 


bending 


moment 


for 


the 


web; 


M 0 : 


Fully 


plastic 


bending 


moment 


for 


the 


flanges; 


M*: 


(u u t v 2 


: ) /4 












N e : 


Fully 


plastic 


membrane 


force 


for 


the 


web; 


N c : 


Fully 


plastic 


membrane 


force 


for 


the 


flanges ; 



GEOMETRIC PARAMETERS 



A: Crushing depth; 

A*: A / H; 

5: Local crushing depth for global bending mode; 

Vertical extent of damage for local denting mode; 
C: C / H; 

a: Angle between first and second hinge lines; 

P: Angle between second and third hinge lines; 
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GEOMETRIC PARAMETERS, CONTINUED 



0 O : Angle of indentation due to crushing; 

0 ± : Angle between regions connected to the i'th hinge 

line; 

1 A : Length of i'th hinge line; 

O: Angle between the first and third hinge lines; 

<J) : Angle between center line and third hinge line; 

y: Angle of depression in global bending mode; 

u f : Stretching of flange; 

u x : Maximum stretching of web; 

e f : Average strain for flange; 

e x : Maximum strain for web; 

e„ g : Average strain for structural unit; 

P B : Crushing force. 
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3. LOCAL DENTING MODE 



3.1. SIMPLIFIED MODEL: Figure 6 shows a photograph of a 
fully clamped plate which was deformed by a crushing load 
applied by a rigid punch to the upper edge of the plate. This 
photograph is from full scale experiments conducted by S. 
Shimizu of Shimushu University in Japan. 1101 Based on this 
photograph and other photographs of damage incurred by tankers 
during grounding accidents, the model for the local deformation 
mechanism shown in Figure 7 was developed. 

The deformation of the web involves both bending and 
stretching. The bending of the web is caused by the formation 
of plastic hinge lines which are shown in Figure 7 . The web 
deforms symmetrically about the point of application of the 
load. Six stationary plastic hinge lines form. Due to 
symmetry, only the three hinge lines on one side of the vertical 
plane of symmetry will be described. The first hinge line forms 
at the junction between the web and the upper flange. The 
second hinge line forms at an angle a below the first hinge 
line. The third hinge line forms at an angle (3 below the second 
hinge line. 

The photographs examined in developing this model showed 
the length of the deformed region in the longitudinal direction 
to be approximately twice the length in the vertical direction. 
This characteristic was included in the model by fixing the 
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angle between the upper edge of the undeformed web and the third 
hinge line at 7i/4 . The vertical length of the deformed zone is 
taken as a parameter C, . 

3.2. DEFORMATION SEQUENCE: Figure 8 depicts the geometry 
of the deformed web. The applied force P B causes the upper edge 
of the web to indent a distance A. The angle of indentation is 
0 O . As A increases the second hinge moves out of plane. The 
angle between the first and third hinges, <E>, decreases. The 
angle between the vertical plane of symmetry and the third hinge 
line, <j>, remains constant at rc/4 . 

The shaded area in Figure 8 illustrates the stretching in 
the deforming web. The attached flange deforms with the web as 
shown in Figure 9. Again, the shaded area illustrates the 
stretching in the deforming flange. 

3.3. LOAD-DEFLECTION RELATIONSHIP: Following the Upper- 
Bound Theorem, the load deflection relationship is derived by 
equating the rate of external work performed in deforming the 
web girder to the rate of internal energy dissipation. In the 
following sections, this derivation is discussed. First the 
rate of internal energy dissipation is derived. Next, the rate 
of external work is formulated. Then, these two rates are 
equated to yield the load-deflection relationship. 
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3.3.1. RATE OF INTERNAL ENERGY DISSIPATION: The 

total rate of energy dissipation consists of two components-- 
that contributed from bending, E 0 , and that contributed from 
stretching, E M . 



E total Zm + E b 



( 2 ) 



3. 3. 1.1. BENDING ENERGY RATE: In general the 

bending energy rate term consists of energy dissipated by the 
continuous deformation field and energy dissipated by the 
discontinuous deformation field. 

E b = f M 0 K dS + t M 0 [6j 1 ± (3) 

J s 1=1 



where S is the surface area of the deforming region; K is the 
rate of change of curvature; 0* is the rate of rotation in the 
i'th hinge line; and 1 A is the length of the i'th hinge line. 

For the local denting model, however, the deformation is 
simplified by modelling all plastic hinges as stationary hinges. 
There is no continuous deformation field. The bending energy 
rate term simplifies to: 

e b - 2 m 0 { rdj l, * [6 2 i j 2 * ie 3 i i 3 ) , (4) 
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where M 0 is the fully plastic bending moment: 




( 5 ) 



The assumption is made that all hinge lengths remain constant 



determined from the geometry of the model shown in Figure 6 . 



The rotation rates, 0 ± , for each hinge line can be related 
to the angle <I> between the first and third hinge lines . For 
each hinge line, a section normal to the hinge line is used to 
determine the rotation rate. 

From the geometry of the model, the following can be 
determined: 



during the deformation process. The hinge lengths can be 



1 1 = c 

1 2 = C / cos (a) 

1 3 = n C 



( 6 ) 



cos 0 X 



cos P - cos a cos $ 
sin a sin $ 



( 7 ) 



cos 0 2 



cos a - cos p cos 
sin P sin 



( 8 ) 
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( 9 ) 



cos 0 3 



- cos a cos P + cos $ 
sin a sin p 



Using the relationship that : 



[ arccos (u) ] 



- u 



SU - u 2 ) 



( 10 ) 



with (7), (8), and (9), the rotation rates can be shown to be: 



ex - 



- 4>(cosa - cosP cos®) 
sin® + c 2 cos® - cos 2 ® 



(ID 



6 , 



® sin® 



^Ci + C 2 cos® - cos 2 ® 



( 12 ) 



0 , 



- ® (cosP - cosa cos® ) 
sin® ^/c x + C 2 cos® - cos 2 ® 



(13) 



where : 



= (sin 2 a - cos 2 p) 
C 2 = (2 cosa cosP) 



(14) 



By noting that the sum of the depression angle, Q e , and ® 
remains constant and equal to k/4, COS (®) , SIN(®), and ® can be 
related to A and A. 
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cosO 



(15) 



= y/2 (C + A) 

2 VC 2 + A 2 



sinO = 

2 VC 2 + A 2 



(16) 



<fc = 



-A C 



(t 2 + A 2 ) 



(17) 



Substituting (11), (12), (13), and (6) into (4) yields: 



E =-2 CM (co3 z g -cosacosPco34>) +/2 (costtcosP-co3 2 gcos<fr) +2sin 2 <I > 

coso sinO + C 2 cos«I> - cos 2 <I> 



(18) 



„ * (C, - C. cosO - cos 2 <I>) 

= -2 C M 0 ^ * — 

cosa sin«I> + C 2 cos$> - cos 2 ® 



(19) 



where 



C 3 = (cos 2 a + v/2cosacosP + 1) 

C 4 = (/2cos 2 a + cosacosP) . (20) 
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3 . 3 . 1 . 2 . MEMBRANE ENERGY RATE : Now consider the 



membrane energy rate term of (2) . Both the attached upper 
plating (or flange) and the web stretch while deforming. 



Therefore : 



# ♦ t 

total ~ flange + web 



( 21 ) 



The rate of membrane energy dissipation is defined as: 

E m = f N 0 t dS = f N 0 k dr\di , (22) 

J s J s 



where e is the velocity strain in the deformed co-ordinates 
(£,T|) and N 0 is the fully plastic membrane force, defined below: 

4 = ' (23) 

dt. 

N 0 = o Q t . (24) 



Integration of (22) yields: 



’lo 

4 = 2 N 0 f u dr] 

0 



(25) 
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where 2u is the velocity with which material points on either 
side of the centerline of the shaded regions shown in Figures 8 
and 9 move apart relative to each other. 



First, consider the stretching involved in deforming the 
flange. Figure 9 shows the trough shape of deformation that the 
flange assumes. Although the stretching of the flange occurs 
across the entire deformed region, the shaded area in Figure 9 
is used to approximate the average strain involved. The amount 
of stretching is constant in the ^-direction, normal to the 
longitudinal girder. 

From the geometry of the model, it follows that: 



72 C 



(26) 



COS0 o 



(cos® + sin®) 



-y/2 ® C (cos® - sin®) 
(1+2 cos® sin®) 



(27) 



Therefore 



b 



flange 




(28) 



o 



= -272 % B { ® 



(1 + 2cos®sin®) 



(cos® - sin®) 
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where N e is the fully plastic membrane force for the flange: 



N o = 



( 29 ) 



Now, consider the stretching involved in deforming the web. 
Again, the shaded region in Figure 8 can be used to approximate 
the average strain for the web. The amount of stretching varies 
along the depth of the web. The displacement u £ at the web 
flange connection and the maximum displacement u x at the 
intersection of the second hinge lines are used to determine the 
average displacement as shown in Figure 10. 



u 



web , avg 



_ i 



/ 



Ctana 



u f + 



(u x - u f ) z 



f 



c 

Ctana 



Ctana 

u x (C - z) 

C ( 1 - tana ) 



dz 



dz 



( 30 ) 



“webavg = \ u x + \ tana u t 



( 31 ) 



# 

^web avg 



1 *1 
— lL + — tana u f 

2 x 2 f 



( 32 ) 
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From the geometry of the model, the following can be determined: 



u. 



^ [ (y/2 coscc - cosP)sin® + cosP cost - cosa] (33) 

[cosp (sin® - cos®) + cosa] 



u. 



= -y/2 ® C 



(cosacosp -cos 2 acos®) 
[cosp (sin® -cos®) + cosa] 



( 34 ) 



Substituting (27) and (33) into (31) and then substituting 
(31) into (25) yields: 



f 

veb 



-n $ c 2 n 0 



(cosacosP - cos 2 acos®) 
[cosp (sin® - cos®) + cosa] 2 

+ sing (cos® - sin®) ^ 
cosa (1 +2 cos® sin®) 



( 35 ) 



3.3.2. RATE OF EXTERNAL WORK: The rate of work done 

in deforming the structural unit is: 

K xc = -Pfl A . (36) 



Rearranging (17) and substituting this into (35) gives: 






~ P B ® 



(C 2 + A 2 ) 

c 



( 37 ) 
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3.3.3. PHYSICAL LIMIT: For the local denting mode, 
a is modelled as a parameter. The choice of a affects not only 
the magnitude of the crushing force for a specified crushing 
depth, but also the range of freedom for the model. The degree 
of freedom for the denting mode is exhausted when the upper 
flange makes contact with the "upper flap" of the web, the 
region between the first and second hinge lines . When this 
occurs, 0! = 7C/2 and : 



A* ^| _ COS 2 « + 2 COsP y/cOS 2 g ~ COS 2 P 

C J limit 2 COS 2 p - COS 2 CC 



(38) 



For approximately a <, n/12, the linear assumption discussed 
in Appendix A is valid. The load-deflection relationship in 
dimensionless terms is: 

Pg = 2 [ (iq - K 2 ) - (K 2 + 2) i|f] 

M o cosa(l - i|r) JUC. 3 - JQ) + (tf 4 - 2) i|r 

+ 2^2 C*#* [ Us - ^ 6 ) - < 39 > 

cos 2 a - 

+ »|i [4 tana + 8 B*t* 2 } 



where the dimensionless parameters are: 



A* = 

B* = 



A 

H ' 
B_ 





(40) 
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and 



A 

C 



( 41 ) 



♦ = - 



The constants K x through K 7 are functions of a and |3 which are 
listed in Appendix A (A6) . 



3.3.4. INITIAL PHASE ENVELOPE: For the local denting 
mode, it is postulated that the longitudinal and vertical extent 
of deformation, is initially zero but spreads during the 
deformation process until the deformed region extends the depth 
of the web (£ = H) . 

Figures 11 and 12 show the local denting solution for 
values of increasing from = 0.1 to = 1.0. The parameters 
for the hull geometry (H* = 100; B* = 75; t* = 1.0) in both 
Figures 11 and 12 are approximate values for a 40,000 DWT USDH 
tanker. Figure 11 depicts the load-deflection curve for a = 
n/12 . The load-deflection curve for a = k /8 is shown in Figure 
12 . 

Both Figures show envelopes for the initial phase of 
deformation. Due to the complexity of the load-deflection 
equation (39) it is difficult to analytically determine this 
envelope. The initial phase, however, may be described 
empirically by a function of the form: 
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(42) 



£b 

Mo 



C 0 (AT 



where the constant C 0 and the exponent "m" depend on both the 
value of alpha and the geometric parameters describing the 
structure . 



3.3.5. STRAIN AND STRAIN HARDENING: From the 
geometry of the model the strain averaged over the entire 
deforming area can be shown to be : 



€ 



avg 



1 









B*t* + tanal e f + 


f 


1 - 

5 




2h' ) 


U h*j 


A 



(43) 

) 



where 



v/2 - cos® - sin® 
(cos® + sin®) 



(44) 



and 



_ (y/2cos« - cosP)sin® + cosPcos® - coscc ( 45 ) 

x [cos(3(sin® - cos®) + cosa] 

Strain hardening effects can be accounted for by adjusting 
the value of a o based on the average plastic strain. For 
example, if the stress-strain relationship for a given material 
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can be described by a power-law relationship similar to the 
Ramberg-Osgood formula C111 : 



o 



o 




(46) 



then, the following refinement of the denting model could 
be made: 



II 


fH 


f ^ 


M ulc 


V € uit J 


\ M 0 



(47) 



where 



% 




(48) 



and P b /Mo is defined by (39) or (A5) . 
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GLOBAL BENDING MODE 



4 . 



4.1. SIMPLIFIED MODEL: This second mode of failure 
involves generation of a three hinge beam-bending mechanism, 
shown in Figure 13. For simplicity, the crushing load is 
applied at mid-span and the length between "hinge points" is 
equal. This mechanism could be adjusted to account for unequal 
lengths between hinge points . 

The "hinge points" are modeled as shown in Figure 14. This 
model is similar to the local indentation model except there is 
no extension of the flange. The vertical extent of deformation, 
£, extends the entire depth of the web/ therefore, C, = H. Also, 
a fourth hinge line forms along the bottom flange at the 
centerline of the "hinge point" mechanism. 

4.2. DEFORMATION SEQUENCE: Figures 13 & 14 show the web 
deforming in the global mode. The applied force P B causes the 
legs of the three hinge mechanism to rotate an angle y. This 
causes the angle <}> for each "hinge point" mechanism to decrease. 
The second hinge line moves out of plane and the angle between 
the first and third hinge lines, O, decreases. This results in 
a local deflection, 5, at the central "hinge point" mechanism. 
This local 5 is part of the total deflection A as shown in 
Figure 13. The shaded area in Figure 15 illustrates the 
stretching of the web involved in deforming the "hinge point" 
mechanism. 
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The range of A* is limited to A* ^ 0.5. For a ship which 
may have a web girder depth on the order of 2 meters. A* = 0.5 
is a significant displacement of the hull. Because the range of 
A* is limited, the extension of the two "legs" of the global 
mechanism is neglected. 

4.3. LOAD-DEFLECTION RELATIONSHIP: The application of P B 
causes the three "hinge point" mechanisms to form. The "hinge 
point" mechanisms are separated by a distance (L/2 - H) . The 
central "hinge point" is designated the primary mechanism. The 
two outer "hinge points" are designated secondary mechanisms. 
The subscripts "p" and "s" are used to designate the primary and 
secondary mechanisms, respectively. As the central mechanism 
deflects, the two "legs" of the overall structure rotate an 
angle y. From Figures 13 and 14, the angles <}> p and <j>. can be 
related to y. 



4>p = 




(49) 



and 



4 >5 = 



TZ_ _ _y 
4 2 



(50) 
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In the following sections, the internal rate of energy 
dissipation for the primary mechanism is formulated as a 
function of y and y, f p (y,y) . In the global bending mode, energy 
is dissipated not only by the primary "hinge point" mechanism 
but also by the two secondary "hinge point" mechanisms. The 
total energy dissipated, therefore, is: 



where f p and f. are the functions formulated in the following 
sections for the total rate of internal energy dissipation for 
the primary and secondary "hinge point" mechanisms, 
respectively. The functions f p and f, are related by (4 9) and 
(50) so that: 



deformation of the web girder in the global bending mode, like 
that in the local denting mode, involves both bending and 
stretching. The expression for the internal rate of energy 
dissipation, therefore, is the same as stated by (2) . 




fp(y.y) + 2f s (y, y) 



(51) 




(52) 



4.3.1. RATE OF INTERNAL ENERGY DISSIPATION: The 
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4. 3. 1.1. BENDING ENERGY RATE: Like the local 

denting model, the bending model consists of six stationary 
hinge lines in the web. In addition, one hinge line forms 
across the bottom flange as shown in Figure 13. There is no 
continuous deformation field. The bending energy rate term, 
therefore, is: 

4p = 2M o { ii ip] + h [& 2 p] + i 3 [6 3p ] } + Mo i 4 [ 0 4p ] (53) 



where M 0 is the fully plastic bending moment for the flange: 



(54) 



The assumption is made that the lengths of the hinge lines 
remain constant. From the geometry of the model, shown in 
Figure 14, the hinge line lengths can be determined: 

h = H 

l z = H/ COSO 

1 3 = s[2 H (55) 

1 4 = H{ 1 - tana) 

1 5 = B 



First, consider the fourth hinge line. The rotation rate 
of the fourth hinge line is simply: 



0 



4 P 



= 2 Y 



(56) 
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Now consider the other hinge lines. The rotation rates, 
[0 A ] , for the first, second and third hinge lines are related to 
the angle O between the first and third hinge lines by the same 
relationships used for the local denting model-- (11) , (12) and 

(13) . This is due to the similar geometry of the two models. 
The relationship of O to A, however, differs from that used for 
the local denting model . 

To relate <E> p to y, the following geometric relationships, 
which are evident from Figure 15, are used: 



p 




(57) 



cos0 o = s[2 sin<J> p 



(58) 



sin0 o = \j cos 2 <|)p - sin 2 <J) p 



(59) 





( 61 ) 
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( 62 ) 



-2y{sin2Yv/sin2Y _ sin 2 2Y + cos2Y ( -j -sin2y) ) 
Vsin2Y~sin 2 2Y v^l - sin2Y + 2sin 2 2Y + 2cos2Y v^sin2Y-sin 2 2Y 



In Appendix B, the E B term is presented as a function of y. 



4. 3. 1.2. MEMBRANE ENERGY RATE: The rate of 

membrane energy rate term is defined by (22) - (25) . For the web, 
the maximum displacement, u x , shown in Figure 15, is given by 
(32) with £ = H. 

The average displacement, u w<ITg for the shaded region in 
Figure 15 is : 



u 



web, avg 




(63) 



Differentiation of (63) and substitution of this into (25) 
yields : 



■“Mr, 



= -s/2 <f> p N 0 H 2 



cosacosp - cos 2 acos4> p 
{ cosP (sin4» - cos® ) + cosa } 



(64) 



4. 3. 1.3. SECONDARY MECHANISMS: The internal 
rate of energy dissipation for the secondary mechanisms follows 
the same function as that for the primary mechanism evaluated at 
y/2 . Both the bending energy rate term and the membrane energy 
rate term for the secondary mechanisms are listed in Appendix B . 
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4.3.2. RATE OF EXTERNAL WORK: The rate of 

external work performed in deforming the web girder is: 

W„ c = P, A . (65) 

The local deflection, 5 , at the central mechanism can be 
described as a function of y by considering the geometry of the 
model . 



6 - yJJP - 2/f z sin 2 <j) J3 = H \/sin2Y (66) 

The total deflection, A, is the sum of the local deflection and 
the deflection due to the rotation of the two legs of the entire 
mechanism. 



Therefore 



A 



A = 6 + - h\ tany 



= 6 + y (-^ - fl) sec 2 y 



(67) 



( 68 ) 



Differentiation of (66) and substitution of the result into 
(68) gives: 




( H (cos 2 y-sin 2 y ) cos 2 y + (L - 2 H) y/cosysiny } 
cos 2 y v'cosysiny 



(69) 
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Substitution of this result into (65) gives: 



W ext = ( -Pfl H y ) Z T 



( 70 ) 



where 



Z r 



A 

Hy 

^ (cos 2 y-sin 2 y) cos 2 y + - 2| Vcosysiny 

2 cos 2 y ycosysiny 



(71) 



4.3.3. PHYSICAL LIMIT: As in the local denting 
mode, a is modelled as a parameter for the global bending mode 
also. The choice of a affects not only the magnitude of the 
crushing force for a specified crushing depth, but also the 
range of freedom for the model. The degree of freedom for the 
bending mode is exhausted when the connection of the web to the 
flange moves out of the Z-Y plane. This occurs when: 



Y llmi C 



4 



arctan 




sin 2 cc H 
sin 2 P ) 



(72) 
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4.3.4. GLOBAL LOAD -DEFLECT I ON RELATIONSHIP : In 

Appendix B, the external rate of work is equated to the total 
rate of internal energy dissipation for both the primary 
mechanism and the two secondary mechanisms . The dimensionless 
equation for P B /M 0 is: 



' B _ 



M „ 



4Z„ 



Cj-QcosG^-cos 2 ®^ 



2y/2H* (cosacosp-cos 2 acosG> s ) 
cosasin^q+qcos^-cos 2 ^ (cosp (sin® p -cos® p ) +cosa) 2 



8 Z e 



C 3 -C 4 cos$ s -cos 2 O s 



+ 2>/2H* (cosacosp -cos 2 acos < I > s ) 
cosasint^Cj+CjCostj-cos 2 ®^ (cosP ( sin$> s -cosG> s ) +cosa) 2 

4 

ZjH* 

(73) 



where 



A* = v/2 (y) 2 + -l] Y 



(74) 



and the a new dimensionless parameter is introduced: 



L * = — 

H 



(75) 



36 



The expression for Z T is given in (71) . The expressions for 
Z p , COS Op, and SIN O p are given by (B2) , (B3) and (B4) 
respectively. The expressions for Z M , COS O,, and SIN O, are 
given by (B6) , (B7) , and (B8) , respectively. The dimensionless 
parameters are listed in (39) and a new dimensionless parameter 
is introduced: 

For small values of y, small angle approximations can be 
employed and higher order terms can be neglected. Because the 
extension of the two "legs" of the global bending mechanism is 
neglected, A* is limited to A* ^ 0.5. Over this range of A*, y 
£ 5 degrees for values of H* and L* which are reasonable for ship 
scantlings . Applying small angle approximations and neglecting 
higher order simplifies the following: 



cos® 



p 




3 

2 



sin® 



p 




(76) 



Z, 




'p 



3 1 



{ y + 2y 2 } 2 
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cos<X> = 


v/2 


3 


2 


sin<I> s = 


& 

o 



1 - (Y> + (Y) 2 + -i (Y) 2 



Z m = 



1 - (Y) 



2(1-y 



(Y) 



( 77 ) 



} 



2 1 
{ Y + 2y 2 - 2y 2 } 2 



Z T 

A* 



y/2 + (L* - 2) (Y) 2 
2 ( Y ) 1 



\/2 ( y ) 2 




1 (Y) 



(78) 



4.3.5. STRAIN HARDENING: Strain hardening effects may be 
included in the global bending model in the same manner in which 
they were included for the local denting model in section 3.3.5. 
For the global bending model, P B /M C in equation (47) is 
represented by (73) . 

4.4. COMBINED BENDING-DENTING MODEL: Figures 16 and 17 
show the combined local denting and global bending load- 
deflection curves for typical parameters for a 40,000 DWT USDH 
tanker. Figure 16 is the load-deflection curve for a = k/12. 
The curve for a = n/Q is shown in Figure 17. Both Figure 
indicate transitions from the local denting mode to the global 
bending mode where the two solutions cross. 
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5 . 



CRUSHING EXPERIMENTS 



5.1. BACKGROUND: Four extruded aluminium double hollow 

sections were crushed for this study. Two experiments were set- 
up to initiate the local denting mode . The other two 
experiments were designed to simulate the global bending mode. 
Appendix C contains descriptions of the experiments . 

The depth-to-thickness ratio for these specimens was 
H/t„ = 20.3, which is at the lower end of the thin-wall member 
spectrum. Also, the material exhibited significant strain- 
hardening properties. Appendix D contains the specification of 
the aluminium specimens that were crushed. 

The extruded aluminium specimens were provided by Reynolds 
Aluminium, Holland B.V. This size of extruded shape was used in 
order to maintain a small scale which could easily be handled 
with standard laboratory machines. The specimens were 300 mm 
long, 127 mm wide, and 50.8 mm deep. The issue of structural 
discontinuity and inhomogeneity due to welds was avoided by 
using extruded forms. 

The double hollow profile was chosen in order to provide 
some restraint against in-plane motion in the transverse 
direction. A multiple cell section such as that shown in Figure 
5 would have been an ideal choice. The double hollow form, 
however, was selected due to availability. 
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5.2. LOCAL DENTING: Figures 18 and 19 show the load- 
deflection curves for the two denting experiments. A photograph 
of the deformed specimen from the first denting experiment is 
shown in Figure 20. 

5.2.1. DESCRIPTION OF DEFORMATION: Due to the significant 
strain hardening effects, it was difficult to see the 
propagation of the deformation during the crushing process. The 
photograph in Figure 20 shows the permanent deformation of the 
first dented specimen with one of the outer webs cut away to 
expose the interior web. The pattern of deformation is very 
similar to that postulated for the simplified model. The angle 
a between the flange-web connection and the first hinge line is 
approximately 15 degrees. The shape of the lower hinge line is 
parabolic rather than straight. The longitudinal length of the 
deformed region, however, is approximately twice the depth of 
the web . 

In both denting experiments, fracture occurred at the 
juncture of the outer two webs and the upper flange . During the 
crushing process, the location of fracture was obscured by the 
rigid punch. Therefore, the crushing depth (A*) at which 
fracture began was not recorded. 

Both denting specimens showed initial signs of fracture 
along the second hinge line. This appears to be a result of the 
fact that a structure cannot assume infinitely small radii of 
curvature. Because the local denting model employs stationary 
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hinges and is formulated for only a discontinuous deformation 
field, this feature is not included in the model. 

The ends of the specimens for both denting experiments were 
clamped to the work surface only at the lower flange. A fully 
clamped condition was not achieved. A small degree of "pull-in" 
in the longitudinal direction was observed during the crushing 
process. The local denting model solution was developed with 
the assumption of fully clamped ends. The amount of energy 
dissipated by stretching in the experiments was probably less 
than that predicted by the simplified model. 



5.2.2. LOAD— DEFLECTION CURVES: The load-deflection curves 
labeled "MODEL" in Figures 18 and 19 were developed from 
Equation (A5) with the parameters listed below. The curves 
shown are the envelopes of the curves produced by (A5) . Each 
specimen was modelled as three structural units with the 
following geometric parameters : 

a = rc/12 H* = 20.3 
B* = 16.9 t* = 1.0 

Although the two outer structural units are channel type 
sections rather than I-beam type sections, it is assumed that 
the two outer webs deform in the same manner as the inner web 
during the crushing process. 
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Strain hardening effects were included as discussed in 
section (3.3.3). The stress-strain curve for the specimens 
which was provided by the manufacturer is included in Appendix 
D. 

As shown in Figures 18 and 19, the local denting model 
over-predicts the energy absorbed in the crushing process. It 
also does not capture the load reduction which occurs at 
approximately A* = 0.10. The model shows an almost linear 
increase in load after the initial phase. Some of the load 
reduction may be due to the lack of fully clamped end 
conditions . 

From the derivations in Appendix A, it can be shown that 
the membrane energy involved in stretching the flange is the 
most significant factor in this region of the load-deflection 
curve. The relative magnitude of the crushing load for the 
model and the general character of the curves, however, show 
good agreement with the experimental results . 



5.3. GLOBAL BENDING MODE: Figures 21 and 22 show the 
load-deflection curves for the two bending experiments. Figure 
23 is a photograph of the deformed specimen from the first 
bending experiment . 
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5.3.1. DESCRIPTION OF DEFORMATION: Again, the 
photograph in Figure 23 shows the permanent deformation viewed 
after experiments were completed. The deformation is similar to 
that postulated for the bending model . The angle a is 
approximately 22 degrees. The parabolic shapes of the hinge 
lines are more pronounced in these specimens than in those of 
the denting experiments . 

In the first bending experiment, fracture occurred at the 
juncture of the two outer webs with the upper flange. The 
second bending experiment was halted at a smaller value of A* . 
The initiation of fracture is evident at the same locations for 
this specimen . 

5.3.2. LOAD-DEFLECTION CURVES: Figures 21 and 22 show the 
load-deflection curves for the bending experiments. The curves 
labeled "BENDING MODEL" and "DENTING MODEL" were developed from 
the models discussed in sections 3 and 4 with strain hardening 
effects included as discussed in sections 3.3.5 and 4.3.5. 

The specimens were modelled as three structural units with 
the following geometric parameters: 

a = 71/ 8 H* = 20.3 
B* = 16.9 t* = 1.0 
L* = 4.4 

The global bending mode was treated as bending of a free 
beam. The contributions from the secondary hinge mechanisms 
were neglected. 
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The theoretical model for the bending mode shows very close 
agreement with the experimental results both in the shape and 
the magnitude of the load-deflection curves. The theoretical 
denting model, however, over-predicts the required force during 
the initial phase prior to transition to overall bending. 
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CONCLUSIONS & RECOMMENDATIONS 



6 . 



The objective of this study was to investigate the crushing 
characteristics of web girders for a USDH structure. The 
crushing experiments that were conducted provide more insight 
into the shape of the deformation that a web girder undergoes 
when subjected to a concentrated lateral collapse load. The 
shapes of the deformed specimens shown in Figures 20 and 23 are 
very similar to the photograph in Figure 6 which shows damage to 
a fully clamped plate subjected to a crushing load. 

The local denting model over-predicted the amount of energy 
absorbed for the denting experiments. A model with a continuous 
deformation field may better predict the initial phase for the 
denting solution. 

From the derivations of the load-deflection for the denting 
model, it can be seen that the energy involved in stretching the 
flange is the most significant term contributing to the total 
rate of internal energy dissipation. The bending energy and the 
energy involved in stretching the web are negligible in 
comparison. Denting experiments with fully clamped end 
conditions would provide better insight into the stretching 
involved in deforming the flange plating. Also, experiments 
using thinner specimens (greater depth-to-thickness ratios) with 
less significant strain hardening characteristics would provide 
more insight into the initial mode of deformation and any 
transition to different modes of deformation. 
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The theoretical load-deflection curve for the global 
bending model closely followed the experimental results for a 
structural unit with free-free end conditions . The shape of the 
load-deflection curve for the global bending model is very 
similar to the shape of experimental results presented by de 
Oliveira 1121 for a horizontally free circular tube subjected to 
a concentrated lateral crushing load. 

The membrane energy involved in the extension of the two 
"legs" of the bending mechanism was neglected for the 
theoretical model. This may be a significant factor for a fully 
clamped structural unit . Crushing experiments for fully clamped 
conditions would be worthwhile. 

The models presented in this paper were based on stationary 
straight hinge lines. Although the actual hinge lines are 
parabolic curves, the models provide reasonable estimates of the 
energy absorbed during the crushing process. For a given set of 
geometric parameters (H*, B*, T*, and L*) , the formulas developed 
for these models could be approximated by simple functions and 
then incorporated in damage prediction models. For example, a 
linearly increasing function could be used to approximate the 
local crushing solution. For the combined denting-bending mode, 
a mean crushing force could be employed. These simple functions 
could be used to describe the springs of the USDH model shown in 
Figure 3. The denting model could be used to describe the 
vertical springs. The global bending model could be used to 
describe the diagonal springs . 
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For this study, several assumptions were made in order to 
keep the analysis of the crushing characteristics relatively 
simple and mathematically tractable. These assumptions may 
require further investigation. 

In particular, the assumption that there is no load- 
interaction between bending moments and membrane forces was made 
to uncouple the energy dissipation terms. The inclusion of the 
interaction surface for the different forces would make the 
problem significantly more difficult, as pointed out by de 
Oliveira 1131 . 

A further refinement of the models presented in this paper 
could be made for the method of including strain hardening 
effects. Strain hardening was accounted for by a relatively 
crude manner which utilized the average strain over the entire 
structure. The strain for each deforming region could be used 
to adjust the flow stress for each region. 

Another factor which requires further investigation is 
fracture. Not only may fracture affect the deformation mode and 
the load carrying capacity of a web girder unit, it may also 
initiate flooding in a ship's double hull. For the models, an 
estimate of initial fracture could be made based using a 
critical strain to rupture criteria for the regions of maximum 
strain . 

For the crushing experiments, fracture first occurred at 
the upper web flange connections for the outer flanges. The 
ultimate strain for the experimental specimens was on the order 
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of 20%. Strain in this connection region should be 
investigated. The extruded shapes had stress concentrations in 
this region due to the geometry shown in Figure 2 6 . A USDH ship 
structure, also, will have stress concentrations due to welded 
joints . 
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APPENDIX A: LOCAL DENTING MODE 



In this section the development and the normalization of 
the load-deflection relationship for the local denting mode is 
detailed. 

A.l. LOAD-DEFLECTION RELATIONSHIP: Equating the external 
rate of work to the internal rate of energy dissipation gives: 



Substituting (37) , (19) , (28) , and (35) into (Al) and making 

dimensionless yields: 




(Al) 



Pb = C! 



2 (C 3 - C 4 cos® - COS 2 ®) 



Mo (C*2 + A* 2 ) 



cosasin® ^/C 1 +C 2 cos®-cos z ® 



+ 4v/2C*tf* [ 



(coscccosp-cos 2 ttcos®) 



{ cosp (sin®-cos®) +cosa } 2 



sing (cos®-sin®) j 
cosa (l+2cos®sin®) 



+ 8v/2 B'C ,Z 



b * c *2 (cos®-sin®) 
(l+2cos®sin® ) 



(A2) 



where 
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cos4> 



= (C * + A*) 

2 



*2 



sin4> = ±± 



s/2 ((« - A*) 

2 v'C* 2 + A* 2 



(A3) 



The dimensionless parameters variables are: 



A* = A ; = Jl . if* = — • 

// H 



S* = — ; t* = 






(A4) 



Substituting (A3) into (A2) yields: 



• s _ 






1 + t 2 



2 [jq (l + t 2 ) (1 + 40 (l+4f 2 ) 1/2 ] 

cosa (l-i|/) ^fK 3 ( 1 +i|/ 2 ) -2i| T+K^ (1+40 (l+4f 2 ) 1/12 



2y/2 C* 

1 + i|r 2 



[K 5 (1+i\> 2 ) -Kg(l+i|r) (l+i|f 2 ) 1/2 ] 



[Kj^-K^ (l+ty 2 ) 1/2 +cosa] 
( (4 C-g. ' tana ^8 iTf 2 ) j ^ (1 + * 2 , 1/2 



(A5) 



*i = 
*2 » 
*3 = 

*4 = 

*5 = 

k 6 = 

*7 = 



2 cos 2 a + 2s/2 cosacosP + 1 
2 cos 2 a + s/2 cosacosP 

1 - 2 (cos 2 a - cos 2 P) 

2s/2 cosacosP 

2 cosacosP 
s/2 cosa 

2 cos 2 p + cos 2 a 



(A6) 
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A * 
C* 



(A7) 



♦ - ' . 



If \(/ < 1, then xj/ 2 « 1 . If this case holds, then Equation 
(A5) can be linearized by neglecting y 2 . Equation (A5) becomes: 



Jj? 

Mo 



+ 



+ 



2 [ Uq - iq) - Uq + 2) x|r ] 

^ cosa (1 - t|r) yj IK 3 -k 4 ) + Uq - 2)H» ; 

' 2/2 C*H* [ (* 5 - K^) W ] > 

; cos 2 a - AT 4 i|r / 

( x|r [4 C*/T tana t8B’t’ 2 ] ) 



(A8) 
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APPENDIX B: GLOBAL BENDING MODE 



B.l. BENDING ENERGY RATE: Substitution of (60), (61), and 

(62) into (11) , (12) , and (13) gives expressions for [0 X ] , [0 2 ] , 

and [0 3 ] in terms of y. Substitution of these expressions and 
(55), and (56) into (53) yields: 



E b = 2 y M 0 H Z p 
+ 



' C 3 -c 4 cos$ X3 -cos 2 $ p > 

^cosasin$ p cos$^cos^ ; 

2 BM c f 



(Bl) 



where C x , and C 2 are given by (14); C 3 ,and C 4 are given by (20); 
and: 




2 {sin2Yv / sin2y-sin 2 2y+cos2Y ( -sin2y) ) 



(B2) 



{sin2y-sin 2 2y} 2 {l-sin2y+2sin 2 2Y+2cos2YY/sin2Y-sin 2 2y} 2 



cosfl> p 



= -^{ (l-sin2y) 



+ (cosy+siny) v/sTn2y } 



(B3) 



sin4>„ - -f{ cos2y - (cosY-sinv )s /sIn^Y } 



(B4) 
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B . 2 . SECONDARY HINGE MECHANISMS: The bending energy rate 

terms and the membrane energy rate terms for the secondary 
mechanisms are equal to those for the primary mechanism 
evaluated at y /2 rather than y. The bending energy rate term 
is : 



E b = 2 y M 0 H Z s 
+ 



/ 



V 



C 3 -c 4 cos$ s -cos 2 $ s 
cosasin® s + cos^~-cos 

B W 0 y 



(B5) 



where 



Z, = 






2 { siny^/siny - sin 2 y + cosy(— - siny) } 

2 

- ^ j 

{siny-sin 2 y } 2 { 1-siny +2sin z y +2cosy\/sinY -sin 2 y } 2 



(B6) 



costf> s = { (1 - siny) + (siny + sin 2 y) 2 } (37) 



sin$ s = { cosy - (siny - sin 2 y) 2 } 



(B8) 
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The membrane energy rate term for the secondary mechanisms 



is : 




-v/2 <t s N 0 H 2 



cosacosP - cos 2 acos$ s 
{ cosP (sin* s - cosQg) + cosa } 2 



(B9) 



B . 3 . LOAD DEFLECTION RELATIONSHIP: The total internal 

rate of energy dissipation consists of the energy rate terms for 
the one primary mechanism and the energy rate terms for the two 
secondary mechanisms. Combining these terms, equating the sum 
to the external rate of work, and making dimensionless yields: 




4Z_ 



Cj-C^COSfcp-COS 2 ^ 



co9asin® pV /C 1 +C 2 cos® p -cos 2 <I> J 



2/2H* (cosacosP-cos 2 acos* s ) 
(cosP (sinfcp-cosflp) +cosa ) 2 



8 zj C 3 -C 4 c os$ s -cos 2 $ s 2</2H* (cosacosP -cos 2 acos* s ) 

z r ( cosa sin* sV /q+C 2 cosa» s -cos 2 4> s (cosp (sin* s -cos* s ) +cosa) 2 

+ 4B*t* 3 

Z^* 

(BIO) 
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where 



A* 



v/2 (Y) ^ 




(Bll) 



The dimensionless parameters are listed in (39) and a new 
dimensionless parameter is introduced: 



L* = 



L 

H 



(B12) 
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APPENDIX C: EXPERIMENT SET-UP 



C.l. DENTING EXPERIMENTS: Figure 24 shows a side-view and 
an end-view of the set-up for the denting experiments. The 
specimens were clamped to the work surface only along the lower 
flange. The work surface (the platen) was raised in a 
displacement control mode and the reaction force P B was measured 
as a function of displacement. A solid steel cylinder 25 mm 
diameter was used as the rigid punch. The solid cylinder was 
used rather than a knife-edge indenter to delay initiation of 
fracture . 

C.2. BENDING EXPERIMENTS: Figure 25 shows a side-view and 
an end-view of the set-up for the bending experiments. A three- 
point bending set-up was used. The specimens were not clamped. 
Instead, bending of a free beam was modelled. The specimens 
were placed on two steel rods of 25 mm diameter that were spaced 
225 mm apart. A third steel rod was connected to the load-cell 
above the specimens. The base was raised in a displacement 
control mode and the reaction force P B was measured as a function 
of displacement . 
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APPENDIX D: EXPERIMENT SPECIMEN SPECIFICATIONS 



D.l. DIMENSIONS: Figure 26 shows a cross-section and the 
dimensions of the prismatic specimens used for the crushing 
experiments . 

D.2. MATERIAL PROPERTIES: The specimens used for the 
crushing experiments were extruded 6061 aluminium which were 
supplied by Reynolds Aluminium Holland, B.V. The material 
properties for the specimens are listed in Table D-l, below. 



TABLE D-l: MATERIAL PROPERTIES 



YOUNG'S MODULUS 


69 


GPa 


YIELD STRESS 


114 


MPa 


ULTIMATE STRESS 


220 


MPa 


ULTIMATE STRAIN 


0.206 


BRINELL HARDNESS 


69 



The material has significant strain hardening properties. 
In order to incorporate this characteristic in the theoretical 
models, the stress-strain curve was modelled as: 



o 




(Dl) 



with n = 0.165; a 0 = 220 MPa; and e ult = 0.206. 
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The points plotted in Figure 27 were obtained from the 



tensile stress-strain curve provided by the supplier, 
shown in Figure 27 is the power-law model which was 
the theoretical model solutions. 



The curve 
used for 
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FIGURES 




Figure 1: The USDH Hull Form 
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Figure 2: Transverse Frame Crushing Model 




Figure 3: Spring Model of USDH 




Figure 4: USDH Section 




Figure 5: Structural Unit 
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Figure 6: Deformation of Clamped Plate 
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t 




Figure 7: Local Denting Model 




Figure 8: Deforming Web Girder 
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Figure 9: Deformed Flange 
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LOCAL DENTING MODEL (ALPHA = PI/12) 



CH* = 100; 8* =75; T* = ID 




Figure 11: Local Denting Mode for a * Jt/12 
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LOCAL DENTING MODEL ( ALPHA = PI/ 8) 



Ct+“ = 100; B* =75; T* = 1} 




Figure 12: Local Denting Mode for a — Jt/8 
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Figure 13: Three Hinge Mechanism 




Figure 14: Hinge Point Mechanism 
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Figure 15: Stretching in Hinge Point Mechanism 
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COMBINED DENTING-BENDING MODEL 



CH*«100; B*«75; T*»1; L*«9; ALPHA=P !/ 12} 




Figure 16: Global Bending Mode for a - rc/12 
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COMBINED DENTING-BENDING MODEL 



CW*=100j B*=75; T*=1j L* = 9j ALPHA=PI/8D 




Figure 17: Global Bending Mode for a « n/Q 
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DENTING #1 



CK*=20. 3; B*=16.9j T*=1J ALPHA = PI/ 123 




Figure 18: First Denting Experiment 
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DENTING #2 



CH*=20.3j B*=16.9j T*=1j ALPHA = P 1/12} 




Figure 19: Second Denting Experiment 
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BENDING #1 



CW*=20.3; B*=16.9j T*«1; L*~4.4j ALPHA=PI/B} 




Figure 21: First Bending Experiment 
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BENDING #2 



C H*=20 . 3j B*-1G.9; T*-1; L*-4 . 4 j ALPHA-PI/B} 




Figure 22: Second Bending Experiment 
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Figure 25: Set-up for Bending Experiments 
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STRESS-STRAIN CURVE FOR ALUMINIUM SPECIMENS 
n * 0.165 



Figure 27 : 




* actua I values 



Stress-Strain Curve for Specimens 
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